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FOREWORD 

This final report represents the second of two volumes of documentation 
for the Generalized Escape System Simulation (GESS) computer program. This 
GESS Programmer”s Manual summarizes program conventions, lists the annotated 
FORTRAN IV program code and describes the procedures necessary to successfully 
use the GESS computer program. Volume I, the GESS User’s Guide, describes the 
elements present and events occurring in typical escape systems, the theory 
and formulation of the simulation model, and the procedures for preparing and 
executing the GESS program and the related ACT and DRAS programs. A portion 
of this work was performed by Ketron, Inc., in accordance with NADC Contract 

& N62269-81-C-0206, Task No. 630-1944. 
| This work was sponsored by the Naval Air Systems Command (AIR<-5312) under 


: Airtask No. A5315312-001D-1531W1100, Work Unit No. WC530. 
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ABSTRACT 


The Generalized Escape System Simulation (GESS) program is a comput- 
erized mathematical model for dynamically simulating the performance of exist- 
ing or developmental aircraft ejection seat systems. The program generates six- 
dimensional trajectory predictions of the aircraft, seat/occupant, occupant 
alone, and seat alone by calculating the forces and moments imposed on these 
elements by the seat catapults, rails, rockets, stabilization, and recovery 
systems included in most escape system configurations. User options are pro- 
vided to simulate the performance of all conventional escape system designs 
under most environmental conditions and aircraft attitudes or trajectories. 
This GESS Programmer's Manual, summarizes program conventions, lists the anno- 
tated FORTRAN IV program code, and represents the second of two volumes of GESS 
documentation. Volume I, the GESS User's Guide, describes the elements present 
and events caren in typical escape systems, the theory and formulation of 


the simulation model, and the procedures necessary to successfully prepare, 


execute, and utilize this and the related ACT and DRAS programs. 
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1.0 INTRODUCTION 


Designing ejection seat systems enabling the safe escape of crew- 
members from high-speed military aircraft is a formidable task. As an adjunct 
to the Navy's on-going escape systems program of system engineering, testing, 
and incident review, a mathematical model has been formulated and developed to 
simulate the operation of any ejection seat-based escape system under most 
realistic conditions from any aircraft. This Programmer's Guide lists the 
annotated code of the Generalized Escape System Simulation (GESS) program, and 
is Volume II of two volumes of GESS documentation. A complete description of 
the history, theory, preparation, and use of GESS is provided in Volume I, 


the GESS User's Guide. 
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2.0 CONVENTIONS 


201 Coordinate Systems 


Trajectories generated by the Generalized Escape System Simulation 
(GESS) program are calculated with respect to specific, “right hand rule,” 
coordinate systems coinciding with center of gravities or other specified ref- 
erence points of the various system elements. Each of the program coordinate 


systems is defined by six (6) degrees of freedom: 


A. Linear: 


1. X =~ forward positive displacement 


nN 
re 
i} 


left positive displacement 


3. Z 


upward positive displacement 
Be Angular: 
4. Yaw (R) - leftward positive rotation 
5. Pitch (Q)- nose downward positive rotation 


6. Roll (P) - right wing downward positive rotation 
These coordinate systems, illustrated in Figure 2-1, are defined as follows: 


e@ Earth-Fixed Coordinate Svstem (EFCS) 

The EFCS is a 3-axis orthogonal coordinate system with origin at 
a fixed point located on or near the surface of the earth. Since all simu- 
lated distances are relatively small compared to the earth's radius, the 
errors associated with neglecting the curvature of the earth's surface are 
considered negligible. The movements and rotations of all simulated system 


elements can be described with respect to (wrt) this fixed coordinate system. 
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Seat Coordinate System (SCS) 


a 


y y Aircraft Coordinate System (ACS) 


Earth-Fixed Coordinate System (EFCS) 


Figure 2-1. Coordinate Systems 
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@ Aircraft Coordinate System (ACS) 


The ACS is a 3-axis orthogonal coordinate system with origin at a 
fixed point, usually the aircraft's center of gravity. All linear and angular 
movements of the aircraft are determined by establishing the position of this 
coordinate system wrt the EFCS. The position of all other system elements, 
prior to and, optionally, after the separation of the seat from the aircraft, 


is determined wrt to the ACS. 


@ Rail Coordinate System (RCS) 


The RCS is a 3-axis orthogonal coordinate system which is defined 
and fixed wrt the ACS. Its origin coincides with the mid-point between the 
rails of the lateral (Y) axis defined by the two lowest rail points in the 
ACS, and its vertical (Z) axis is parallel to the rails. The rail system 
serves to restrict the motion of the seat to a specific path as it moves out 


of the aircraft cockpit. 


@ Seat Coordinate System (SCS) 


The SCS is a 3-axis orthogonal coordinate system with origin usually 
corresponding to the initial location of the RCS origin wrt the seat. The 
center of gravity (0G) of either the seat/occupant (S/0) combination or the 
seat/alone (S/A) is referenced wrt the SCS, allowing the trajectories of these 
system elements to be tracked wrt either the ACS or the EFCS. Points of appli~ 
cation of the various forces acting on the S/O or S/A are also referenced wrt 


the SCS, allowing the calculation of moments and rotation around the respective 


CGs. 
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® Occupant Alone Coordinate System (OACS) 

The OACS is a 3-axis orthogonal coordinate system with origin at 
the center of gravity of the seat occupant after the occupant separates from 
the seat. The OACS is currently defined only by the three linear degrees of 
freedom because the disjointed non-rigidity of the occupant alone makes any 
angular determination difficult, if not irrelevant. Provisions have been made 
for the future incorporation of OACS angular degrees of freedom if such 


information is necessary at a later date. 


e@ Thrust Vector Control Coordinate System (TVCCS) 


The TVCCS is a 3-axis orthogonal coordinate system with origin at 
the common intersection of the rocket line of thrust and the centerlines of its 
gimbals. This coordinate system establishes the orientation of the rocket 
line of thrust, which is modified dynamically as part of the vertical-seeking 


maneuver. 


Coordinate systems and vector directions should be carefully con- 


sidered when preparing simulation inputs. All output reports are referenced 


accordingly. 


2.2 Variable Identification 
Several variable naming conventions are used to identify variables 
used in the program. Refer to the program annotation for specific variable 


descriptions. These naming conventions are summarized below. 


SA - seat alone 

OA - occupant alone 
so - seat/occupant 
NPTS - number of points 
RK - rocket 
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WGHT - weight 

x - X axis 

> 4 - Y axis 

Z - Z axis 

P - roll rotation 

Q - pitch rotation 

R - yaw rotation 

VEL - velocity 

PORO - porosity 

REC - recovery parachute 
DRO/DR - drogue parachute 
DRT - DART 

POS - position 

IGN ~ ignition 

CAT - catapult 

TVC - thrust vector control 
CG ~- center of gravity 


The variable names used in the program have been selected to facili- 
tate the programmer's understanding. There is limited correlation between these 


variable names and the variable names used in the program formulation. 


2.3 Integration Arrays 


The large number of equations and saved values associated with many 
of the trajectory integrations in GESS resulted in the development of a highly 
structured integration array logic. This logic maintains the current value of 
each integration equation, in addition to all intermediate or previous equation 


values for both the Runge-Kutta and predictor-corrector methods. 
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Thus, the initial four time steps can be integrated by the Runge-Kutta routine, 


and subsequent 


integrations 


performed by the predictor-corrector routine. 


The first value of each array indicates the number of equations in- 


volved in the integration. 


The remaining array values are then ordered, 


accordance with the number of equations, as described in Table 2-1. 


Array Location 
1 


2 (14+N) 
(2+N) (1+2N) 
(2+2N) (1+3N) 
(2+3N) (1+4N) 
(2+4N) (1+5N) 
(2+5N) (1+6N) 
(2+6N) (1+7N) 
(2+7N) (1+8N) 
(2+8N) (1+9N) 
(2+9N) --(1+10N) 
(2+10N)--(1+1L1N) 
(2+11N)--(1+12N) 
(2+12N)-~(1+13N) 
(2+13N)--(1+14N) 
(2+14N)--(1+15N) 


(2+15N)--(1+16N) 


Integration Array Values 

Description 
number of equations (N) 
results of integrations - equation values 
equation derivatives 
current Runge-Kutta K-values 
intermediate Runge-Kutta K-value summations 
intermediate equation values 
equation values at (time-At) 
equation values at (time-2dt) 
equation values at (time-34t) 
equation values at (time-44t) 
equation derivatives at (time-At) 
equation derivatives at (time-2At) 
equaticn derivatives at (time-3At) 
predictor values at currect time 
predictor values at (timetét) 
corrector values at current tine 


corrector values at (time+at) 


in 
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2.3.1 Aircraft Trajectory Array. The aircraft trajectory array, TRAJAC, 


is based upon the !2 primary trajectory variables described in Table 2-2. 


Table 2-2. Aircraft Trajectory Array Description 


Description 


number of aircraft trajectory equations 
= 12 if tracking aircraft 
= 0 otherwise 


Array Location 
TRAJAC(1) 


TRAJAC( 2) position of aircraft in EFCS 7 


TRAJAC( 3) y position of aircraft in EFCS 


TRAJAC(4) z position of aircraft in EFCS 


TRAJAC(5) velocity of aircraft in ACS 


TRAJAC(6) velocity of aircraft in ACS 


TRAJAC(7) velocity of aircraft in ACS 


TRAJAC(8) (roll) angular position of aircraft ft 


TRAJAC(9) 


(pitch) angular positon of aircraft ft 


TRAJAC( 10) 


(yaw) angular position of aircraft ft 


TRAJAC( 11) (roll) angular velocity of aircraft in ACS 


TRAJAC( 12) 


(pitch) angular velocity of aircraft in ACS 


TRAJAC( 13) (yaw) angular velocity of aircraft in ACS 


TRAJAC(14) derivatives and other equation values associated 
through with the above equations as detailed in Table 2-1 
TRAJAC( 193) ‘ and with N=12. 


t - integration result; not used for describing angular position. 
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2.3.2 Occupant Alone Trajectory Array. The occupant alone trajectory 


array, TRAJOA, is based upon the 12 primary trajectory variables described 


in Table 2-3. 


Table 2-3. Occupant Alone Trajectory Array Description 


| Array Location Description 
7 
? TRAJOA( 1) number of occupant alone trajectory equations 


= 0 before S/O separaton 
= 12 after S/O separation 


TRAJOA( 2) x position of occupant alone in EFCS 

TRAJOA( 3) y position of occupant alone in EFCS 

TRAJOA(4) z position of occupant alone in EFCS 

TRAJOA( 5) x velocity of occupant alone in ACS 

TRAJOA(6) y velocity of occupant alone in ACS 

TRAJOA(7) z velocity of occupant alone in ACS 

TRAJOA(8)* p (roll) angular position of occupant alone 
TRAJOA(9)* q (pitch) angular positon of occupant alone 
TRAJOA( 10)* r (yaw) angular position of occupant alone 
TRAJOA(11)* p (roll) velocity of occupant alone in OACS 
TRAJOAC 12)* q (pitch) velocity of occupant alone in OACS 
TRAJOA(13)* r (yaw) velocity of occupant alone in OACS 
TRAJOA( 14) derivatives and other equation values associated 
through with the above equations as detailed in Table 2-1 
TRAJOA( 193) and with N=12. 


* ~ All angular positions and velocities and their derivatives are currently | 
permanently set = 0; array space has been maintained to facilitate future 


incorporation of O/A angular tracking. | 
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2.3.3 Seat Alone Trajectory Array. The seat alone trajectory array, 


TRAJSA, is based upon the 12 primary trajectory variables described in 


Teble 2-4. 


Table 2-4, 


Array Location 
TRAJSA( 1) 


TRAJSA( 2) 
TRAJSA( 3) 
TRAJSA(4) 
TRAJSA( 5) 
TRAJSA(6) 
TRAJSA( 7) 
TRAJSA(8) 
TRAJSA(9) 
TRAJSA( 10) 
TRAJSA( 11) 
TRAJSA( 12) 
TRAJSA( 13) 
TRAJSA( 14) 


through 
TRAJSA( 193) 


Seat Alone Trajectory Array Description 


Description 
number of seat alone trajectory equations 
= 0 before S/O separaton 
= 12 after S/O separation 
x position of seat alone in EFCS 
y position of seat alone in EFCS 
z position of seat alone in EFCS 
x velocity of seat alone in SCS 
y velocity of seat alone in SCS 
z velocity of seat alone in SCS 
p (roll) angular position of seat alone t 
q (pitch) angular positon of seat alone t 
r (yaw) angular position of seat alone ft 
p (roll) angular velocity of seat alone in SCS 
q (pitch) angular velocity of seat alone in SCS 
t (yaw) angular velocity of seat alone in SCS 
derivatives and other equation values associated 


with the above equations as detailed in Table 2-1 
and with N#12, 


t - integration result, not used for describing angular position. 
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2.3.4 Seat/Occupant Trajectory Array. The seat/occupant trajectory array, 


TRAJSO, is based upon the 12 primary trajectory variables described in Table 


2-5. 
Table 2-5. Seat/QOccupant Trajectory Array Description 

Array Location Description 
TRAJSO( 1) number of seat/occupant trajectory equations 

= 12 before S/O separaton 

= 0 after S/O separation 
TRAJSO( 2) x position of seat/occupant in EFCS 
TRAJSO( 3) y position of seat/occupant in EFCS 
TRAJSO(4) z position of seat/occupant in EFCS 
TRAJSO(5) x velocity of seat/occupant in SCS 
TRAJSO(6) y velocity of seat/occupant in SCS 
TRAJSO(7) z velocity of seat/occupant in SCS 
TRAJSO(8) p (roll) angular position of seat/occupant f 
TRAJSO(9) q (pitch) angular positon of seat/occupant ft 
TRAJSO( 10) r (yaw) angular position of seat/occupant f 
TRAJSO(11) p (roll) angular velocity of seat/occupant in SCS 
TRAJSO( 12) q (pitch) angular velocity of seat/occupant in SCS 
TRAJSO( 13) rt (yaw) angular velocity of seat/occupant in SCS 
TRAJSO( 14) derivatives and other equation values associated 
through with the above equations as detailed in Table 2-l 
TRAJSO( 193) and with N#12, 


t = integration result, not used for describing angular position. 
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2.3.5 Parachute Trajectory Array. The parachute array, TRAJCH, is based 
upon 6 primary trajectory variables for each of the 3 parachutes incorporated 


into the GESS model, as described in Table 2-6. 


Table 2-6. Parachute Trajectory Array Description 


Array Location Description 
TRAJCH(1,1*) number of parachute i* trajectory equations 


= 6 when applicable 
= 0 otherwise 


TRAJCH(2,1*) x position of parachute i* in EFCS 

TRAJCH(3,1i*) y position of parachute i* in EFCS 

TRAJCH(4,1*) z position of parachute i* in EFCS 

TRAJCH(5,i*) x velocity of parachute i* in EFCS 

TRAJCH(6,i*) y velocity of parachute i* in EFCS 

TRAJCH(7 ,i*) z velocity of parachute i* in EFCS 

TRAJCH(8,i*) derivatives and other equation values associated 
through with the above equations as detailed in Table 2-1 
TRAJCH(97,1*) and with N = 6, 


* -i = 1 for first drogue parachute or slug/container 
= 2 for second drogue parachute 


= 3 for recovery parachute 
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2.3.6 Thrust Vector Control Array. The thrust vector control array, 


TVCEQS, is based upon the 14 primary control variables described in Table 2-7. 


Table 2-7, Thrust Vector Control Array Description 


Array Location Description 
TVCEQS(1) number of thrust vector control (TVC) equations 


= 14 during TVC 
= 0 otherwise 


TVCEQS(2) x\* 

TVCEQS(3) x2* 

TVCEQS(4) x3* 

TVCEQS(5) x4* 

TVCEQS(6) x5* 

TVCEQS(7) x6* 

TVCEQS(8) x7* 

TVCEQS(9) xg* 

TVCEQS( 10) xg* 

TVCEQS(11) *10* 

TVCEQS( 12) 11% 

TVCEQS( 13) *12* 

TVCEQS( 14) x13 

TVCEQS(15) ; x14" 

TVCEQS(16) derivatives and other equation values associated 
through with the above equations as detailed in Table 2-1 
TVCEQS(255) and with N = 14, 


* - for details see User's Manual, Section 3.5.3, and Subroutine VERTSK, 


page 5-303. Note that xy = TVCVALS(1i) and xy 2 TVCDERV(i). 


| 
E 
| 
| 
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2.3.7 Quaternion Arrays. The 4 quaternion arrays, QUATAC, QUATOA, QUATSA, 


and QUATSO, are based upon the 4 primary quaternion array variables described 


in Table 2-8, 


Array Location 
QUATxx*(1) 


QUATxx* (2) 
QUATxx*(3) 
QUATxx*(4) 
QUATxx* (5) 
QUATxx*(6) 


through 
QUATxx* (65) 


AC for 
OA for 
SA for 
SO for 


* - xx 


Table 2-8. Quaternion Array Descriptions 


Description 
number of xx* quaternion equations 
= 4 ,£ xx*® trajectory is being tracked 
= 0 otherwise 
NOxx (See User's Manual, Eq. 3.a) 
Atxx (See User's Manual, Eq. 3.b) 
A2xx (See User's Manual, Eq. 3.c) 
A3xx (See User's Manual, Eq. 3.d) 
derivatives and other equations associated with 


the above equations as detailed in Table 2-1 and 
with N = 4, 


aircraft trajectory 
occupant alone trajectory 
seat alone trajectory 
seat/occupant traiectory 
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3.0 AERODYNAMIC COEFFICIENT TABLES 


During the execution of the GESS program, equally spaced aerodynamic 
coefficients are required for computing the aerodynamic forces acting upon the 
seat/occupant system. These coefficients are usually derived through wind 
tunnel experiments using the subject ejection seat(1)* the forces and moments 
data produced by these tests are used as inputs for the Aerodynamic Coefficient 
Table (ACT) program(2) , which creates the necessary coefficient tables as func~ 
tions of the seat orientation and velocity. The ACT program provides an effic- 
fent method for creating, modifying, and storing these aerodynamic coefficient 
tables on random access files for retrieval by GESS during simulation execution. 


The ACT program was closely modeled after the RFWTHR program (3), 


The annotated code listings of the ACT program are provided in Appen- 
dix A. These listings have been extracted from Reference 2 because of the im- 
portance of having proper aerodynamic coefficient tables for the successful ex- 
ecution of GESS. Detailed instructions for using the ACT program are given in 


the Volume I GESS User's Guide. 


* All references may be found in Section 6.0. 
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4.0 EXECUTION JOB STREAM 


The GESS program is executed by submitting to the computer a series 
of operating system commands known as the “job stream”. The job stream can be 
entered line-by-line interactively. However, it is usually more convenient to 
prepare a job stream, or “submit file", in advance of execution and “batch” 
submit this file to the computer for execution. The interactive commands to 
batch submit a job stream are: 

GET, filename 


SUBMIT, filename 


Another print site can be specified if the main computer center printer is not 


desired. 


SUBMIT filename ,EIl=print site identifier 


If no printed output is desired: 


SUBMIT, filename,N 


Figure 4-1 presents a recommended job stream for GESS execution using 


the CDC KRONOS operating system. This job stream provides for: 


. the loading of the prepared input file (GESSI) 


. the loading of the ACT program-generated aerodynamic coefficient 
tables (AERO4) 

. the loading and execution of the pre-compiled program binary code 
(GESSB) 


. the optional creation or replacement of program~-generated plot- 


ting data files (GESST and GESSZ) 


NADC-8 4068-60 
the creation of a single, inclusive, indirect access output file 
containing a mirror copy of the original input file and all 
requested program-generated reports (GESSO) 
- A “day file” record of the submitted job execution (GESSDAY) 
. An error recovery logic that saves the day file and any prior out- 


put, in the eventuality of an execution error 


The recommended job stream is intended for execution runs generating 
a nominal number of output reports and/or simulating a relatively few seconds 
of simulated trajectory. Should a large number of reports and/or a long simula- 
tion period be desired, the output (and, possibly, the plotting files) should be 
created either as a direct access file, or as several indirect files. Refer to 


Reference 4 for details on direct file creation on the KRONOS system. 
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/ JOB 
GESS(CB200000,T75) 
ACCOUNT ( XXXXXX, YYYYYY) 
ASSIGN,MS, OUTPUT 
GET(TAPE1L=KGESL) 
GET (TAPE2AERO4) 
GET(GESSB) 
MAP (OFF) 
GESSB. 
REPLACE (TAPE42*KGEST ) 
REPLACE ( TAPE41*KGESZ) 
GOTO,1l. 

EXIT. 

1,REWIND,«. 
SKIPF(DAYFILE) 
COPYSBF(TAPE1, OUTPUT) 
COPYBF (TAPES, OUTPUT) 
COPYBF(TAPE6, OUTPUT) 
COPYBF(TAPE7 , OUTPUT) 


. 


COPYBF (TAPE4O, OUTPUT) 
REWIND, OUTPUT. 
REPLACE (OUTPUT =KGESO) | 
GOTO, 2. 

EXIT. 

2,DAYFILE,GESSDAY. 

REPLACE, GESSDAY. 
EXIT. | 
/EOR 
/EOL 


Figure 4-1. GESS Execution Job Stream for CDC KRONOS/OS 
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5.0 GENERALIZED ESCAPE SYSTEM SIMULATION 


(GESS ) PROGRAM 
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ANNOTATED CODE LISTINGS 
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